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Abstract. We report on our recent work about the description of a meson gas below the chiral phase 
transition within the framework of Chiral Perturbation Theory. As an alternative to the standard 
treatment, we present a calculation of the quark condensate which combines the virial expansion and 
the meson-meson scattering data. We have also calculated the full one-loop elastic pion scattering 
^vq [ amplitude at finite temperature and we have unitarized the amplitude using the Inverse Amplitude 

Method in order to reproduce the temperature effects on the mass and width of the p and a 
resonances. Our results show a clear increase of the thermal p width, as expected from previous 
analysis. The results for the a are consistent with Chiral Symmetry Restoration. We comment on 
the relevance of our results within the context of Relativistic Heavy Ion Collisions. 

o 

^ INTRODUCTION AND MOTIVATION 

The recent development of the Relativistic Heavy Ion Collision (RHIC) program is 
one of the main motivations to study hadronic matter under extreme conditions of 
temperature T and density. Here we will consider the meson gas formed when the 
plasma created after one such RHIC has expanded and hadronized, being in a state 
where the chiral symmetry is restored. There are strong indications that one could 
observe the medium effects on such a system. For instance, the dilepton spectrum shows 
an anomalous behaviour for invariant masses near the p mass flj]]. The flatness of the 
spectrum is consistent with a modification of the mass and width of the p's which 
have time to decay inside the plasma, so that their spectral function acquires thermal 
corrections due to the interaction with the hot and dense hadron gas [[![ |3], |J |5], [| [7], |[ p . 

In such a system, external momenta and temperature are small compared to the chiral 
symmetry breaking scale ~ 1 GeV. The relevant degrees of freedom are then the 
lightest mesons and the interactions among them are best described by a low-energy 
QCD effective theory based on chiral symmetry. The most general framework com- 
prising the QCD Chiral Symmetry Breaking pattern SU L (N f ) x SU R (N f ) -> SU v (N f ) is 
Chiral Perturbation Theory (ChPT) [(T^, [TTJ, [12|, [T3]] where any observable can be calcu- 



X 



lated as an expansion in p/A%, p denoting a meson mass, momenta or the temperature. 
Thus, ChPT provides model-independent predictions, just by fixing a few low-energy 
constants (LEC). This program has included the calculation of thermodynamic observ- 



ables such as the free energy density and the quark condensate, as we will briefly review 
below. Throughout this paper, we will neglect finite baryon density effects, as ideally in 
the central rapidity region formed after a RHIC. 

Since ChPT is intended to provide a systematic low-energy and low-T perturbative 
expansion, we do not expect it to reproduce resonances. This is strongly related to 
the fact that ChPT only satisfies unitarity in a perturbative fashion. Over the last few 
years, there has been a lot of work devoted to enlarge the ChPT applicability range by 
using unitarization methods, i.e, imposing exact unitarity requirements, like the Inverse 
Amplitude Method (IAM) [[Rp or approaches based on Lippmann-Schwinger or Bethe- 
Salpeter equations [|T5j] . These methods provide a good agreement with the experimental 
phase shifts and they are able to generate resonances, like the p and the o for the 
SU (2) chiral symmetry. In addition, they can be extended to include coupled channels, 
describing successfully all the meson-meson data and resonances in the SU (3) case, up 
to 1.2GeV[[rj],[l7|]. 

What we will show below is that only requiring chiral symmetry and unitarity one 
can also describe successfully the thermal behaviour of the p and a resonances. For 
that purpose, one needs first to calculate the thermal pion scattering in ChPT, which 
has been done in [|18|]. We shall discuss the main features of such thermal amplitude 
below. Then, by using the IAM extended to finite T, one can construct a nonperturbative 
thermal unitarized amplitude which, in particular, describes the behaviour of the thermal 
p and o [|19|]. We will present the results for the thermal mass and width of the p and o 
in that approach, as well as for the T-dependence of the effective g pK7t vertex. The main 
implications of our results in the context of RHIC and chiral symmetry restoration will 
be also discussed below. 



THE MESON GAS AND CHIRAL SYMMETRY AT FINITE T 

WITHIN CHPT 

For the reasons explained above, it is important to provide an accurate description 
of the low-T meson gas in thermal equilibrium. For instance, the signature of chiral 
symmetry restoration at T c ~ 150-200 MeV should be observed in the thermal evolution 
of the order parameter, the quark condensate (qq)(T) from below the transition point. 
As we have just discussed, ChPT provides a model independent description of the 
thermodynamic observables, based only on chiral symmetry. The only extra ingredient 
is the temperature, which is treated as an &{p) parameter. 

The first calculations of the pion gas within ChPT go back to [KJ, where (qq)(T) 
and the thermal dependence of f n (T) were calculated up to &{T^ (one loop) in the 
chiral limit. That result already showed a behaviour compatible with chiral symmetry 
restoration. In [|T]j a thorough analysis up to <ff(T 6 ) was performed, including the free 
energy, the quark condensate and an estimate of the thermal effects of free kaons and 
etas. The &(T ) corrections to f K (T) have been analyzed in [ |2~2"1 ] where it has been 



shown that beyond <ff(T 2 ) one has to consider separately the space and time components 
of the axial current, so that there are two independent which, in addition, can be 
complex. In fact, their imaginary part is proportional to the in-medium pion damping 



rate while their real parts are related to the deviations of the pion velocity from the 
speed of light. Other analysis of the thermal pion dispersion law can be found in [E3L 1231] . 



The analysis of typical nonequilibrium effects such as explosive pion production after 
a RHIC can be also studied within the ChPT context p5|]. It should be pointed out 



that many of these properties have also been investigated using specific models for low- 
energy QCD. The most popular is the 0(A) model, which reproduces a critical behaviour 
already in the mean field approach. Apart from introducing the o explicitly, conventional 
perturbation theory in the 0(A) model breaks down, which has been dealt with at finite 
T using different nonperturbative approaches [ |2"6"| , |Z7| ]. 

When calculating thermodynamic quantities such as the pressure or the quark con- 
densate from ChPT, the usual approach is to use the Feynman rules of Thermal Field 



Theory Q28Q to the order considered. This is particularly cumbersome in the three fla- 
vor case if one wishes to include the full dependence on temperature, quark masses and 
SU (3) interactions. An alternative |29| , |30| ] is to perform a virial expansion of the 



pressure as [31] 



/3P = £^.(7X.+£ U 7 (r)^+i^.(r)^.j (i) 

where i = n,K,r]. This is a dilute gas expansion in the fugacities £ ; = exp(— mjT). The 
binary interactions between the different species show up in [£TJ, ^|, [50p: 

B\f = ^k- dEE%(E/T) £ (2/+ 1)(27+ l)S^ JS {E), (2) 

where K { is the first modified Bessel function and 8y are the elastic scattering ij — > ij 

phase shifts at T=0 (chosen so that 8 = at threshold) of a state ij with isospin /, 
angular momentum J and strangeness S. What makes the virial expansion useful is that 
the T dependence on thermodynamical observables up to T ~ 200 — 250 MeV can be 
obtained just from the T = phase shifts, which have been calculated for all possible 
meson-meson interactions in SU (3) ChPT to one loop. They can be found for instance in 
[[T7[]. Note that for the pressure one could even use the experimental phase shifts directly. 
However, the quark condensate is given by the derivative of the pressure with respect to 
the quark mass and therefore the analytic dependence of the 5' J with the different meson 
masses is needed. 

For temperatures much below 150 MeV, massive states like kaons and eta are ther- 
mally suppressed, typically by the Boltzmann factors exp(— M K /T) and, in principle, 
the suppression is even stronger for the interactions among them and with pions, as 
(P shows. For low T it is then reasonable to treat those states as free particles, as it 
was done in [[Hp. Moreover, at low T the integrals in ([Z]) are dominated by the phase 
shifts at threshold. However, for higher temperatures, the effect of massive states be- 
comes increasingly important and, furthermore, the dependence of the interactions with 
the pion mass can be large so that their contribution to the quark condensate becomes 
sizable [BDfl . When the effect of the strange states is taken into account, there are two 



main questions that can be analyzed: on the one hand, the effect on the non-strange con- 
densate (uu) = (dd) — (qq) /2 (in the isospin limit) of adding another flavor. Since the 
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FIGURE 1. The thermal evolution of the (qq) and (ss) condensates in ChPT 



number of degrees of freedom increases, so it does the entropy and one expects that the 
collective state is closer to "disorder". This would imply a decrease of the critical tem- 
perature, as it is indeed observed in lattice calculations [f32"l], although with gluons and 
massless quarks. On the other hand, one can study the thermal evolution of the strange 
condensate (ss). Since m s ^> m one expects that the thermal evolution for (ss)(T) 
is slower than for the non-strange condensate. Note that the quark masses play a simi- 
lar role as external magnetic fields in ferromagnets, so that increasing the external field 
means that it takes more thermal energy to restore the symmetry. Let us remark that, in 
contrast to the lattice, the physical masses are easily incorporated in our approach. 

We report here on the virial expansion calculation in SU (3) ChPT that has been done 
in [[301] . The main results for the quark condensates are summarized in Figure |T} First, 
it should be pointed out that the curves are plotted at most up to the point where the 
condensates vanish. The ChPT condensates do not vanish above those "critical" values, 
since they have been obtained from a perturbative series in the temperature. In fact, the 
results should be trusted only in the region showed in the graphs, as explained above. The 
thermal evolution of (qq)(T) is shown in different approximations, namely, considering 
only pions in SU (2) or SU (3) (their tiny difference comes from the 0(p 4 ) phase shifts), 
adding free kaons and etas [|]J and, finally, considering the full SU (3) interactions 
[30|]. Note also that these interactions (basically only TtK and tzt] are important at these 
temperatures) yield a larger effect than naively expected. The reason is that they depend 
strongly on m K , which is more sensitive to m ud than m K or m^. Taking into account 
the numerical values of the LEC's with their errors, one gets a reduction of the melting 
temperature of T^ su ^ - 7^)^(3) = 21+* 4 MeV, in remarkable agreement with the 
chiral limit lattice calculations [[3"2|], taking into account that we have used the actual 
meson masses. In addition, we estimate T^ su< ^ = 21 1+^ 9 MeV []30[]. The second 
effect, is also clearly seen in Figure [I]: The thermal evolution of the strange condensate 
from the broken phase is much slower than the non-strange one. From Figure HJ, we see 
that there is still about 80 % left for (ss)(T) / (ss)(Q) when (qq) (T) / (qq) (0) has already 



melted. Finally, one may wonder about the effect of calculating the integrals in @ with 
ChPT, whose applicability does not extend to infinity. Indeed, when (g) is evaluated 
[ [301 ] with the phase shifts unitarized with the coupled channel IAM [ |T7| ] , which have a 
much larger applicability range, the numerical results only change very slightly. Thus, 
the main conclusions remain the same, since, as we have already commented, the main 
contribution to the integrals comes from the low-energy region and the IAM agrees with 
ChPT at low energy, improving only the high energy behavior. 



PION SCATTERING AT T ^ IN ONE LOOP CHPT 

If the pion gas is dilute enough, i.e., at very low temperatures, it is reasonable to 
assume that the only dependence of pion scattering with T can be accounted for through 
the Bose-Einstein distribution functions and one can ignore the T-dependence of the 
scattering amplitudes. However, that dependence could be important in several contexts. 
For instance, it has been suggested that an enhancement of pion scattering in the scalar 



channel near threshold could be an indication of chiral symmetry restoration Q27|]. In 



addition, a previous calculation [H, 131] of thermal nn scattering in the Nambu-Jona- 



Lasinio model, shows a singular behaviour of the scattering lengths at some critical 
temperature, which may be related to a Mott transition. Furthermore, if one wishes to 
extend to finite T the fruitful unitarization program in ChPT to describe resonances, 
the full calculation of the %% scattering amplitude to one loop and the extension of 
perturbative unitarity are essential ingredients. 

For the above reasons, it is important to provide a model-independent description of 
pion scattering for temperatures T well below the chiral scale A%. This can be achieved 
in ChPT. The calculation of the thermal scattering lengths in ChPT to one loop has been 
done in whereas the calculation of the full thermal amplitude in one loop ChPT has 



been recently carried out [ |18| ] . We will report here the main results and features of that 
work. 

There are two formal aspects regarding the calculation of the scattering amplitude at 
finite T that are worth pointing out. The first one is that we are considering the thermal 
amplitude defined by taking T = asymptotic pion sates, the T dependence coming from 
the four-point function, calculated using the Feynman rules of Thermal Field Theory 



in the imaginary-time formalism [EBp . Then, one can perform an analytic continuation 



from discrete frequencies co n = 2nnT to real energies E as i(o n — ► E + ie. Such analytic 
continuation corresponds to the retarded four-point function in the so called real-time 
formalism [[35|]. The retarded functions have a causal and analytic structure [ |36| ] suitable 
to extend perturbative unitarity at T ^ 0. The same definition of thermal amplitude has 
been used in [1,0]. 

The second point is that the loss of Lorentz covariance inherent to the thermal formal- 
ism (due to the choice of the thermal bath rest frame) means that any two-body scattering 
amplitude with four-momenta k { k 2 — * k^k 4 will depend separately on the variables S , 

|S|, T , |T|, U and |U|, where S = k i + k 2 , T = k l — k 3 and U = k x — k 4 . At T = 0, 
the amplitude depends only on the Mandelstam variables s = S 2 ,t = T 2 , u = U 2 and 
any %% scattering amplitude can be related to that of n + n~ — > n°n , called A(s,t,u), 
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FIGURE 2. One loop diagrams contributing to Kit scattering 

by isospin and crossing transformations. At T 7^ and since the temperature does not 
modify the interaction vertices, the crossing symmetry still holds, but now in terms of 
S, T and U. Therefore, any tzk thermal amplitude can be written in terms of the thermal 
A(S, T, U; j8). The loop diagrams are the same as for T = and are given in Figure 0. 
The final result for the thermal amplitude can be written as: 



A(S,T,U;/3)=A 2 (S,T,U)+Af / (S,T,U)+A^(S,T,U;/3)+Af(S,T,U;/3), (3) 

where A 2 is tree level contribution coming from the &{p ) lagrangian (the nonlinear 
sigma model) while A 1 ™ 1 contains both the tree level 0(p 4 ) contributions plus poly- 
nomials coming from the renormalization of the loop integrals. Both A 2 and A p ° l are 

temperature independent. The A^ m term represents those contributions from diagrams 
a,b,c in Figure ^[ They yield the correct analytic structure and will ensure perturbative 
unitarity in all channels. Finally, the contribution A4 accounts for tadpoles like those 
in diagrams d,e in Figure □ plus terms coming from diagrams a,b,c proportional to the 
tadpole integral. 

The detailed results for the different contributions above can be found in |Jl~5|] and we 



do not give them here for brevity. As a first check of consistency, we recover the T — > + 



limit of [JTTp. Furthermore, when the thermal amplitude is projected into partial waves 
djj of definite isospin / and angular momentum J (defined in the center of mass frame 
where the pions are at rest with the thermal bath) we also agree with the results given in 



||34J] for the scattering lengths. 

Another important check of consistency, which will be crucial for our analysis in the 
next section, concerns the imaginary part of the partial waves and perturbative unitarity. 
At T = 0, unitarity constraints the partial waves, for s > 4m| and below other inelastic 
thresholds, to satisfy 

lma(s) = o(s)\a(s)\ 2 , (4) 
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FIGURE 3. Interpretation of thermal bath contributions to the pion scattering phase space. In the 
process A, the medium enhances the two-pion production, while process B represents the absorption 
of the initial pions by the thermal bath. 



where o(s) = a/1 — 4-rn^/s is the two-pion phase space, whereas the ChPT series only 
satisfies unitarity perturbatively, i.e., 



lma 2 (s) 



0, 



lma 4 (s) 



cr(s) \a 2 (s) | 



(5) 



It is possible to generalize the relation @ to the case of any one-loop elastic scattering 
amplitude at finite T . The derivation is given in [ jTE| ] and follows closely the analysis in 
[ |37| ] of the discontinuity in the self-energy of a particle decaying in the thermal bath. 
The result for the thermal perturbative unitarity relation is: 



I I L 

Ima 2 (5)=0, Ima 4 (5;/3) = <7 r (S ) \a 2 (s) | , S > 2/7%, 



where 



o T {E) 



a(E 2 ) 



1 + 



exp(j6|£|/2)-l 



(6) 



(7) 



and it has been assumed that only %% states are available in the thermal bath. Thus, o T 
is nothing but the thermal phase space, whose origin can be physically interpreted by 
writing 1 + 2n B (E/2) = [1 +n B (E/2)} 2 - n 2 B (E/2) where n B (x) = (exp(x/T) - 1)) _1 
is the Bose-Einstein distribution function. Written in this way, the first term represents 
the enhancement of phase space due to the presence of pions in the medium, while the 
second term accounts for the absorption of the two initial pions by pions in the bath. 
This is schematically depicted in Figure [31. We have checked explicitly that the relation 
(Eh holds for our thermal %n scattering amplitude. Moreover, the partial waves at T ^ 
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FIGURE 4. Temperature evolution of the phase shifts 8 U for IJ = 1 1 , 00, 20. 



can be analytically continued to the s complex plane [ [191 ] and they display the same 
analytic structure as the T = one, i.e, they have a right unitarity cut in the real axis 
above s > 4m| (the discontinuity across the cut is given by @) and a left one for s < 0. 
We remark that our results remain valid when the density of states with more than two 
pions is small. This is equivalent to neglecting higher powers of n B and is therefore 
similar to the dilute gas approach. For small energies and temperatures, that is, strictly 
within ChPT, the dilute gas approach is consistent, as we have also seen in the previous 
section. When the range of energies is extended, as we will do in the next section, one 
should bear in mind that the range of T where our approach is valid is such that the 
Bose-Einstein factors still remain small. 

Finally, we have plotted the thermal phase shifts in Figure 0. These results deserve 
some comments. First, we observe that the absolute value of the phase shifts in all 
channels increases with T, while their sign (i.e, the attractive or repulsive nature of 
each channel) is preserved. Recall that the phase shifts are related to the real part of 
the amplitude as 8 T ~ o T (y/s) [a 2 (s) + Rea 4 (s, j8)l . The dominant contribution to the 
phase shifts thermal enhancement is given by the thermal phase space factor a T , while 
the T-dependence of the real part of the amplitude is rather weak at low T. In particular, 
we do not observe any significant thermal enhancement for the real part of a 00 , which 
would be interpreted as chiral symmetry restoration p7[]. However, it must be stressed 
that we are just considering the perturbative amplitude, valid only for low T. In the next 
section, we will consider a nonperturbative extension of the amplitude, which does show 
a behaviour compatible with chiral symmetry restoration in the 00 channel. 

Another important comment regarding our results shown in Figure J| concerns the 1 1 
channel, i.e, the p channel. Recall that, following the hypothesis of resonance saturation, 
increasing 8 n is equivalent to increasing the ratio T p f^/M 5 p [[IT], ^8)]. Therefore, our 
results are consistent with a thermal increase of the rho width coming mostly from 
thermal phase space and an almost constant rho mass, at very low T 0, @], §||. This 
behaviour will be confirmed by our analysis in the next section, where we will find also 
important corrections at higher temperatures. 



THE THERMAL p AND o 



The purpose of this section is to show that only from chiral symmetry and unitarity one 
finds a thermal evolution of the masses and widths of the p and o at rest consistent with 
previous analysis and with the dilepton spectrum observed in RHIC. Unlike the models 
where the resonances are included as explicit degrees of freedom, we will start from the 
model-independent pion scattering amplitude in one loop ChPT calculated in the previ- 
ous section and, imposing exact unitarity, we will construct a nonperturbative amplitude 
whose poles in the second Riemann sheet in the 00 and 1 1 channels correspond to the o 
and p respectively. 

The exact unitarity relation (Q) implies that any partial wave satisfies a = 1/ (Rea -1 — 
ia) on the real axis below inelastic thresholds. A unitarization method is just one way 
of approximating Re a 1 . The IAM uses the one-loop ChPT and thus ensures that exact 
unitarity is exactly satisfied and, at the same time, the low-energy predictions of ChPT 
are preserved. The IAM unitarized amplitude for one channel reads a IAM = a|/ (a 2 — a 4 ) 
[[Bp and coincides formally with the [1,1] Pade approximant in the f^ 2 expansion. 

We have already sketched in the introduction the virtues of the IAM formula at T = 
for the description of resonances and the data for higher energies. At T ^ 0, we have 
seen in the previous section that, perturbatively in ChPT and for a dilute gas, the partial 
waves satisfy thermal unitarity @-(|7|). Therefore, following the same steps as for the 
T = case and motivated by the success of the IAM approach, we will consider the 
unitarized IAM thermal amplitude: 

a lAM {s-T)= a ' {s) (8) 
a 2 (s) -a 4 (s;T) 



which satisfies the exact elastic unitarity condition 



lma IAM (s-T)=a T (s 



a IAM (s;T] 



(9) 



The first hint that (|8|) provides a proper description of thermal resonances comes from 
the following simple argument. The behaviour of the 1 1 partial wave in the real axis near 
s = Mj, the p mass squared, can be described by a Breit-Wigner parametrization (valid 
for narrow resonances Y T <C M T ): 

a BW (s;T)= *J {S 1 - (10) 

where R T (s) can be related to the pnn effective vertex (see below). Now, compare (|8|) 
with ( |T0| ) near s = Mj. First, one gets Rea 4 (My.) = a 2 (Mj), which defines the resonance 
mass and, second, from the unitarity relation @, we have T T M T = —R t {Mj)o t (Mj). 
Therefore, if the thermal corrections to R T and to M T were much smaller than those to 
r r (R T ~ R Q and M T ~ M Q ) we would simply get 

r r ^r [i+2 % (M /2)], (ii) 




which is the behaviour expected at very low T for a p at res t [|4| , |3} [8|] that we had 
already anticipated in the previous section. As discussed in [|37|7@] and also in the 
previous section, this result accounts for the difference between the direct decay p — > Tin 
and the inverse one %% — > p which is allowed in the thermal bath and is responsible 
for the dilepton production. Note that in order to derive < JTT| ) we have neglected the 
T-dependence in Rea 4 (5;r). Therefore, by using the complete result for the thermal 
one-loop amplitude discussed in the previous section, we will be able to find, for 
higher temperatures, the T dependence on M T , R T , possible deviations from the low- 
T behaviour (]TT|), and even more importantly, the thermal evolution of both the p and a 
poles in the complex plane, which is the consistent way to generate resonances within 
ChPT. Recall that a Breit-Wigner description is not valid for the a. According to our 
previous discussion, the upper limit in T to which our approach is valid will be dictated 
by the condition n B (M/2, T) < 1 where M is the mass of the resonance described (p or 
a). This gives roughly T < 300 MeV for the p and T < 180 MeV for the a. 

In the above discussion, it is crucial to obtain the analytic continuation of the thermal 
amplitude to the complex s-plane. The details can be found in [|19|]. The result shows 
the same right and left cuts structure as the T = amplitude, the discontinuities across 
the cuts being T-dependent. Once such analytic continuation has been obtained, it can 
be continued to the second Riemann sheet, using (g), as a u (s;T) = a IAM {s;T)/[\ — 
2ia T (s)a IAM (s;T)]. 

Let us first show the results for 8 n for different temperatures, depicted in Figure [5|. 
The SU(2) LEC (see [ JTT| ] for their definition) we have used are = —0.3, l 2 = 5.6, 
Z 3 = 3.4 and 7 4 = 4.3 and are obtained by fitting the T = scattering data, which yields 
M Q = 770 MeV and r o = 159 MeV. We see clearly the broadening of the p as T 
increases. This is confirmed by the evolution of the thermal poles, which is shown in 
Figure The p width increases with T while the p mass decreases slightly. The a pole 
deserves some comments. We see that for low T the o width increases and its mass 
decreases slightly, following a similar behaviour as the p, i.e, mostly dictated by the 
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FIGURE 6. Evolution of the a and p poles with the temperature. 



thermal space factor o T . However, for T > 125 MeV, M G decreases more rapidly and 
T a starts to decrease. This behaviour has an interesting explanation in terms of chiral 
symmetry restoration. On the one hand, since the T-dependence of m K (T) is rather weak 
up to T ~ 200 MeV [ [241 , Wh , the decrease of M a (T) points towards M a — ► m n and 
that decreasing is stronger as T increases, unlike the p mass. On the other hand, as 
T increases and M a approaches 2m % from above, both the direct o — > 2tt and inverse 
111 — > o decays become kinematically disallowed, so that is reduced. At low T, 
the decrease of the mass is much weaker and the phase space contribution dominates, 
making r,j grow. Similar results and discussion for the a can be found in [ |27[ [39]]. 

Finally, in Figure [7| we have collected our results for the p. First, we clearly observe 
a significant deviation (for T > 100 MeV) between the full result (obtained either from 
the pole position in Figure |5| or from the phase shift in Figure |5J) and the naive thermal 
phase space prediction (fTT]), stressing the importance of having a full T-dependent ChPT 
description of %% scattering. Moreover, although M p changes little, consistently with 
Vector Meson Dominance [Q, |]], our results show a sizable slight decrease of the p 
mass for T > 150 MeV, which seems to be favored by phenomenological analysis of 
RHIC dilepton data In Figure [7] we have also plotted the effective pun vertex, 

defined from R T in fllOD as R T = gj (4m|— Mj) /487T, from the VMD pnn coupling 
[0. 0] with a thermal p (g Q ~ 6.2). At low temperatures g T <g (g^/gQ — 0.9991) in 
agreement with the chiral analysis in []7|]. For higher temperatures g T grows, although 
the thermal corrections are much smaller than those at finite baryon density p0[]. 
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CONCLUSIONS 

We have reviewed some basic aspects about the ChPT description of the meson gas 
formed after a Relativistic Heavy ion Collision. In particular, we have shown that one 
can combine the virial expansion with ChPT in order to obtain the thermal evolution 
of the quark condensate below the chiral phase transition in the dilute gas regime. That 
approach allows to use directly all the available information on meson-meson scattering 
for three light flavors in ChPT. Contrary to lattice studies, in this approach physical 
masses are easily accounted for. The main effects of considering the strange sector in the 
condensate are, first, that the strange condensate melts much more slowly than the non- 
strange one and, second, that the melting temperature for the non- strange condensate is 
sizably smaller with three flavors than with two, and that, surprisingly, the kK and KT] 
interactions provide a largest contribution to this effect than free kaons or etas. 

We have also shown the results of a recent calculation of the nn scattering amplitude 
at finite T in one loop ChPT. The partial waves satisfy an extended version of perturba- 
tive unitarity, where the only change with respect to the T = case is that the phase space 
is thermally enhanced by the presence of pions in the thermal bath. Thus, the imaginary 
part of the thermal amplitude has a neat physical interpretation in terms of absorption 
and induced emission of pion pairs in the thermal bath. The thermal phase shifts are also 
enhanced with T, keeping their attractive or repulsive nature, mainly dominated by the 
phase space factor. 

The thermal pion scattering amplitude can be unitarized using the IAM, following 
similar steps as for T = 0. The unitarized partial waves can then be analytically contin- 
ued to the complex energy plane and their poles in the I = J = 1 and 7 = 7 = channels 
correspond to the p and the a. This approach provides a description of the thermal ef- 
fects for a p and a at rest, only from chiral symmetry and unitarity. We have found that 
T p increases significantly with T. For low T, M p and the effective vertex g pnK remain 
almost constant, which is consistent with neglecting the T-dependence in the real part 




of the amplitudes and considering only the contribution from the thermal phase space. 
However, for higher T, the full T-dependence in the amplitude has to be taken into ac- 
count. Thus, for T > 100 MeV, T p acquires significant corrections from the pure thermal 
space approximation and g pnK increases, while for T>150 MeV, M p shows a sizable de- 
creasing behaviour. As for the o pole, although its interpretation as a resonance is much 
less clear than the p, the thermal behaviour of both M a and T , as obtained from its 
associated pole, can be understood from chiral symmetry restoration. 

Our results agree with several theoretical analysis and are consistent with phenomeno- 
logical studies of RHIC dilepton data. We have only used chiral symmetry and unitarity 
as guiding principles, without including the resonances as explicit degrees of freedom. 
It would be interesting to apply our chiral unitary approach to the production of dilep- 
tons from thermal %% annihilation near the p energy, in order to be able to provide more 
accurate predictions about the observed dilepton spectrum in RHIC. This is just but one 
of the possible directions that we will pursue in the near future. 
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